Kinetic roughening in two-phase fluid flow through a random Hele-Shaw cefl 
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A nonlocal interface equation is derived for two-phase fluid flow, with arbitrary wettability and 
viscosity contrast c = (jUi — ^2)/(p.i + M2), in a model porous medium defined as a Hele-Shaw cell 
with random gap &o -I- 5h. Fluctuations of both capillary and viscous pressure are explicitly related 
to the microscopic quenched disorder, yielding conserved, non-conserved and power-law correlated 
noise terms. Two length scales are identified that control the possible scaling regimes and which 
scale with capillary number as £\ ~ 6o(cCa)~^''^ and ^2 ~ feoCa"^. Exponents for forced fluid 
invasion are obtained from numerical simulation and compared with recent experiments. 
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The displacement of a fluid by another in a porous 
medium is a problem of fundamental interest in nonequi- 
librium physics as a paradigm of interface dynamics in 
disordered media jE], ||. Experiments on bead packs in 
Hele-Shaw cells [^in particular, have stimulated con- 
siderable theoretical efforts, but the problem has con- 
sistently revealed itself rather elusive |^. More re- 
cently a new surge of interest has arisen with the recog- 
nition of the inherently non-local character of the prob- 
lem as a key ingredient |Q, and the realization of a 
new series of experiments in Hele-Shaw cells with ran- 
dom gap [H 1^, 1^. Roughening exponents of the pro- 
posed nonlocal equations have been explored by means 
of Flory-type scaling arguments [Q and phase field sim- 
ulations [|[ ^ . While the specific properties of noise are 
known to be crucial to determine the universal aspects of 
interface roughening, fluctuations are usually modeled at 
a phenomenological level, and including only local capil- 
lary effects. Noise related to the non-Laplacian viscous 
pressure due to quenched disorder in the permeability 
has been so far neglected. While this may be justified 
for imbibition experiments other situations, such as 
forced fluid invasion, do require a quantitative assessment 
of this point. In addition, it would be desirable to have 
a unified formulation for general conditions of viscosity 
contrast c = (/xi— ^2)/(mi+M2) and wettability given the 
rich variety of phenomena that the experimental evidence 
has unveiled as a function of those parameters |^ . 

Here we address the general problem of fluid displace- 
ment in a Hele-Shaw cell with random gap, as a sim- 
ple model of a porous medium. This model system has 
the great advantage that no coarse-graining procedure 
must be invoked in the theoretical formulation, thus al- 
lowing us to derive ah initio a general and complete in- 
terface equation, quantitatively accurate, with explicit 
dependence on 'bare' parameters, and including all noise 
sources. On the experimental side, the system is also 
appealing since a direct control of the disorder is locally 
possible on the microscopic scale Q . 

A complete description of interface fluctuations must 
contain three basic physical effects of a porous matrix on 



its motion, namely local variations of (i) capillary pres- 
sure, (ii) permeability, and (iii) available volume. Dif- 
ferent but not independent noise terms must thus be 
generated through distinct physical mechanisms from the 
unique source of randomness. Fluctuating capillary pres- 
sure can be directly related to gap variations Sb = b — ba 
where 60 is the mean value, in terms of the Young-Laplace 
condition for the pressure jump across the interface 
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where k is the curvature in the cell plane, a is surface 
tension and 9 is the contact angle, cos 9 — 1 meaning 
perfect wetting of the (invading) fluid 1. The effect of 
viscous pressure fluctuations, however, is far less obvious 
due to the inherently nonlocal character of the interplay 
of (ii) and (iii) in the response of the fluid flow to gap 
fluctuations. We base our analysis on the assumption 
that, for sufficiently smooth gap variation (i.e. \Vb\ ^ 1) 
Darcy's law for a Hele-Shaw cell H holds locally as 



[bo + Sb{x,y)f 



Vp. 
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In a capillary tube of lateral size h at fixed injection 
pressure, from Darcy's law a relative velocity fluctua- 
tion scales as 5v/v ~ 25b /b (larger permeability implies 
less resistance to flow) while at fixed flow injection it is 
exactly the opposite, 5v/v ~ ^2db/b (mass conservation 
slows down the flow if there is more volume available). 
In an actual disordered medium the solution of the whole 
pressure field will thus be required to sort out the effective 
flow conditions at each location. A direct consequence of 
(iii) is that the 2d effective flow, must give rise in gen- 
eral to a non-conserved interface equation, precisely to 
account for volume conservation in 3d. In our case, 3d 
incompressibility implies that the 2d flow will satisfy 

V • (6v) = 0, (3) 

where the gap acts effectively as a variable density. 
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Eqs. (g) and (g) imply that the pressure field is non- 
Laplacian. In order to obtain a closed bulk equation 
which we can project onto the interface degrees of free- 
dom we treat the pressure perturbatively in |V5|. We 
thus split the pressure field as p = po + Sp and keep only 
the order |V6| in Sp to obtain 



\7^Sp^ 



3V6 



V^Po = 
• Vpo = 0, 



(4a) 
(4b) 



where we have neglected higher orders consistently with 
the fact that they have also been omitted in assuming 
local Darcy flow. The lowest order Eq. (^^ can be solved 
as the unperturbed problem j^, with the modified 
boundary condition Eq. (|l|) which contains all capillary 
effects, while the simplest boundary condition Sp — is 
then required for the Poisson Equation (^) . In terms of 
the Laplacian Green's function Sp then satisfies 



dsG{x - x{s),y- y{s)) 



dSp 
dn 



J dx'dy'G{x ~x',y ~y') • Vpo 



(5) 



The free-boundary problem is thus defined by Eq. (g) 
specified at the interface, with p — po + Sp, and the 
boundary conditions at infinity. Here we focus on the 
case of forced fluid invasion, where a fixed velocity V is 
imposed at infinity and fix > fj.2- We introduce the di- 
mensionless quenched noise as — b'^[l + ^{x, y)]. Noise 
originated respectively from 6b in Eq. (^, from Sb in 
Eq. (^) and from Sp in Eq. (^) will be called respec- 
tively capillary, permeability and bulk noises. 

Concerning the scaling properties of the interface, we 
are interested in the lowest order approximation on the 
interface deviation from planarity, which is relevant in a 
Renormalization Group (RG) sense. Our result for the 
interface equation in Fourier space for two-fluid displace- 
ment under constant injection velocity V takes the form 



1 dh{k) 
V dt 



= S{k) - c\k\ [1 + ihkf] h{k) + Nh{k) 



^-{l+i2\k\)^h[k) + nLR{k,t) (6) 

where the lengths £i and (.2 are defined in terms of 
the capillary number Ca = 12(/ii -I- ^2)V/a as li — 
6o(cCa)~-^/^ and £2 — 26o cos6'Ca^^, and where Mh{k) 
denotes the leading (quadratic) nonlinearities 

/oo 
dq[l - S{kq)\\q\[l + {liqf]Hk - q)h{q) 
-00 

(7) 

where S is the sign function. ^h{k) is the Fourier trans- 
form of ^{x, h{x)) and the term n^ulk, t) is a long-ranged 



correlated noise of the form 
1 



^LR{k,t) = 



Ml + M2 



^iin-{k,t) + ^j,2^+{k,t) (8) 



with 

n^{k,t) = \k\ J dxdy^{x,y + Vt)e-'^''^y\^\Q{±y), (9) 

where Q is the step function. 

Note that the long-ranged term ClLR{k, t) enters effec- 
tively as an annealed (explicitly time-dependent) noise 
(see discussion below). Eqs. (||-||) constitute our central 
result. Note also that in this formulation we have as- 
sumed weak noise so that multiplicative noise terms of 
order h£_ or nonlinear in ^ have been neglected 

The linear deterministic part of Eq. (^ is well 
known The complete set of deterministic nonlin- 

earities can be obtained systematically using the weakly 
nonlinear expansion developed in Ref. [|l^. These in- 
clude the familiar local terms considered in Ref. 0| , such 
as (V/i)^, but also nonlocal terms in real space. 

The capillary fiuctuations give rise to the conserved 
(area-preserving) noise term proportional to |fc|C/i- This 
contribution is associated to the second term in Eq. (|l|) 
and is generated exactly as the usual capillary term 
— \k\^h{k) which comes from the first term in Eq. (j^) 

. A conserved noise term of this form was phenomeno- 
logically argued in Ref. [||. On the other hand, the non- 
conserved noise term proportional to results from the 
trivial lowest order contribution of permeability noise of 
the form VS_h{k), plus a nontrivial local term of opposite 
sign coming from the expression of the bulk noise defined 
as 5vi_ ~ ~j2il^in '^hich takes the form 



Sv^{k) 



W 
"2" 



{-Uk) 



+ ClLRik,t) 



(10) 



We now proceed to sketch the derivation of the bulk 
noise. For simplicity we will consider the one-sided case 
{c ~ 1). Our derivation for c ^ 1 follows the formulation 
of Ref. but is more involved and will be presented 
elsewhere ||ll[. Neglecting orders ^dxh, Eq. (0) reads 



da;' In [{x - x' f + {h{x) - h(x')f\ Sv^(x') = 
dx' {In [{x - x'f + (/i(.t) - h{x')f\ i(x', h(x')) 



^^K-')-y') U,.A (U) 



where the integral on y of the rhs of Eq. (^ has been 
integrated by parts and it has been assumed that the 
noise vanishes at infinity, S^(x,y — > —oo) — 0. It has 
also been applied that ^Vpo - ~|j|f and that the 
Green function has the form Q,(x,y) = — (47r)~-'^ ln[a;^ -t- 
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y^]. With the change of integration y = y' — Vt, and to to the mean interface position |Q, Eq. ( [Til ) then reads 
lowest order in the interface deviation h — Vt with respect 



dx' In \x — x'\5v(_{x') 



dx' \hi\x - x'\^(x' ,h) + I dy 



[x — x'Y + y- 



(12) 
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Eq. dl^ ) can be exphcitly solved for 5v(^ using distribution 
Fourier calculus to yield the result Eq. ([lO|). 

The quenched noise will be typically correlated on a 
microscopic scale a. For ka <^\, effectively white. If 
{£.{^:y)£.{^' ,y')) = A(5(x — x')S{y — y') then we have 

{nLnik,t)nLRik'A')) - A|fc|5(fc + fc')e-l'=l^l*-*'l (13) 

SO CIlr scales as \k\^^^ and introduces long-range mem- 
ory. Accordingly, \ow-k behavior is dominated by the 
local part of bulk noise. This implies that 3d conser- 
vation overcomes permeability at low-fc, giving rise to an 
overall non-conserved noise with the same sign as the cap- 
illary noise. Although direct computation of both terms 
in the bulk noise shows that the local part is typically 
larger, it is unclear to what extent neglecting the long- 
range term may miss important details of local interface 
pinning which may eventually affect the scaling. Fur- 
thermore, for ka 1, both local and nonlocal parts of 
the bulk noise are comparable but then the annealing 
approximation may not be justified, and a more careful 
analysis of bulk noise based in Eq. ( pi] ) may be neces- 
sary |11 . The local and non-local contributions to bulk 
noise may be of the same order in other situations. For 
instance, persistent noise ^ = ^(x) yields Sv^ = 0, with 
non-conserved and conserved noises opposing each other. 

One of the salient features of Eq. (^) is that the prob- 
lem has two characteristic length scales. £i controls 
the well-known crossover between (deterministic) capil- 
lary and viscous forces [^. The second length £2 is a 
newly identified one which defines a crossover between 
conserved and non-conserved noise. For general viscosi- 
ties and wetting conditions the two length scales are ar- 
bitrary and define a variety of possible scaling regimes 
and crossovers depending on their relative size. Experi- 
ments of a wetting fluid invading an inviscid one (c = 1 , 
cos 9 = 1) typically have £2 ^ £i- Concerning the nonlin- 
earities, power counting arguments show that quadratic 
terms are relevant in the RG sense only if the viscous 
term c\k\ is absent, that is, for \k\£i ^ 1. This capillary- 
dominated regime will always be observed at large k and 
short times provided £2 ^ £i- A crossover to the viscous- 
dominated regime will eventually occur for \k\ii ~ 1 and 
then the quadratic nonlinearities do become irrelevant. 
A second crossover within the viscous regime, from con- 
served to non-conserved noise, will occur for \k\£2 ^ 1- 



The explicit knowledge of the bare coefficients of the 
nonlinear terms is helpful to assess the validity of the 
linear approximation. Numerical simulation of typical 
cases show that the effect of nonlinearities in the capil- 
lary regime is not appreciable in reasonable simulation 
times. The capillary regime, which is relevant to experi- 
ments [p[, is thus well described by 



1 dh{k) 
V dt 



{hkfh{k) + ]-£2ih{k) 



(14) 



~ 1 /2 ■ 

with £1 = 60 Ca ' . Since the quadratic nonlinearities 
vanish for c = (for symmetry), we may refer to the 
growth described by Eq. ( |T^ ) , as the universality class of 
'symmetric fluid invasion' (SFI) [ p^ . 

We will now study the SFI scaling by numerical sim- 
ulation of Eq. (|lj) and also the subsequent crossovers. 
The quantities of interest concerning the scaling prop- 
erties are the root mean square of the interface height 
fluctuations, W, and the structure factor S{k,t). The 
width grows with time as W{t) ^ before saturation, 
and the saturation width scales with system size L as 
Wsat ^ L"^- The roughness exponent a can also be ob- 
tained from the relation S{k,t) ~ fc"-'^"^" for long times. 
For further reference, we provide the scaling exponents 
of Eq. (|l|) for two cases that are exactly solvable [ pT| , 
namely, for persistent noise S, = C i^) 1 ^^^d annealed noise 
^ = ^{x,t) which, assuming (5-correlations, give respec- 
tively a = 3/2, /? = 1/2 and a = 0, /3 = 0. 

Typical results for Eq. (|lj) are shown in Figs. |l| and 
|. We choose V = I, a = 0.0625, £1 = 50 and £2 = 3000, 
which satisfy the criterion a ^ £1 <C ^2- For these values, 
we have obtained a roughness exponent a = 1.2 ± 0.05 
(very close to the value reported in Ref. (|] ) and a growth 
exponent f3 = 0.68 ± 0.02. The scaling of the correlation 
function G{1, t) has been found to be fully consistent with 
Q^ioc = 1, as expected from the superrough value of a > 
1. Note that the scaling of the power spectrum S{k,t) 
at large k corresponds to the case of persistent noise, 
a = 3/2: short segments of the interface 'feel' effectively 
a persistent noise for long enough time intervals. 

An increase (by an order of magnitude) in the value of 
V with the subsequent variation of £1 and £2 modifies the 
scaling behavior: for low values of k the scaling of S{k, t) 
yields a = 0, the value obtained for annealed noise, while 
for larger k the observed scaling is essentially the same 
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as in Fig. |[ Hence, the effectively noise acting on the 
interface at low-A; for large V is annealed. The value of 
the exponent /3 is observed to decrease with increasing 
V, consistently with the trend towards an effective an- 
nealed noise, for which (3 = 0. On the other hand, if V is 
decreased by an order of magnitude, the scaling observed 
in Figs. |l| and || is essentially unchanged. 

The values of Ca and V we have used are of simi- 
lar magnitude to the ones of the experimental work of 
Ref. |8|, and the exponents reported there for large in- 
jection rates, a ~ for small k and a ~ 1.3 for large k 
are fully consistent with our results for large V. How- 
ever, this is not the case for small values of V. This 
discrepancy could be attributed to non-darcyan effects 
associated to the sharp edges in the gap variations of the 
experiment, which may modify pinning properties. 

After the first crossover to the viscous regime with 
conserved noise, we find that the interface is not rough, 
Wsat{L — > oo) = const. After the second crossover to 
non-conserved noise, we get a = P = 0, with W growing 
logarithmically both with time t and system size L. 

As a final remark let us mention that the case of spon- 
taneous imbibition (constant pressure conditions) can be 
described with our formalism with an appropriate time- 
dependence V{t) ~ except close to pinning, when 
the zero-mode fluctuations must be carefully worked out. 

We expect that the physical insights and the predictive 
power of the theoretical framework here presented may 
be useful to reinterpret data and design new experiments, 
in particular in the yet unexplored range of parameters. 
The extent to which it may be applicable to more realistic 
porous media, with appropriate coarse-grained parame- 
ters, however, is an open question which deserves further 
study. 
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Terms of this kind coming from capillary and permeabil- 
ity can be included systematically |l^, but those from 
the bulk cannot be obtained explicitly. 
Unlike local noise terms, the 'annealed' approximation is 
justified in Sv^ as a leading order, because of its integral 
form. This means that the noise acting on a region where 
the interface is pinned does vary because the interface is 
moving elsewhere and it couples through the bulk. 
With respect to RG flow of c this is a saddle fixed point. 
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FIG. 1: Interface width W as a, function of time, for system 
sizes L = 32, 64, 128, 256, 512 and 1024. The straight line is 
a fit with a slope P = 0.68. 
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FIG. 2: Structure factor for a system with L = 256. The data 
are for t = 0.5 (lower curve) to t = 12.0, and time interval 
At = 0.5. The straight line with slope —3.3 {a = 1.2) is a fit, 
and the other straight line has slope —4 (a = 3/2). 



